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ABSTRACT 



GRB 090417B was an unusually long burst with a T90 duration of at least 2130 s 
and a multi-peaked light curve at energies of 15-150 keV. It was optically dark and 
has been associated with a bright star-forming galaxy at a redshift of 0.345 that is 
broadly similar to the Milky Way. This is one of the few cases where a host galaxy 
has been clearly identified for a dark gamma-ray burst and thus an ideal candidate 



5_i ■ for studying the origin of dark bursts. We find that the dark nature of GRB 090417B 
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cannot be explained by high redshift, incomplete observations, or unusual physics in 
the production of the afterglow. Assuming the standard relativistic fireball model for 
the afterglow we find that the optical flux is at least 2.5 mag fainter than predicted 
by the X-ray flux. The Swift/XRT X-ray data are consistent with the afterglow being 
obscured by a dense, localized sheet of dust approximately 30-80 pc from the burst 
along the line of sight. Our results suggest that this dust sheet imparts an extinction 
of Ay > 12 mag, which is sufficient to explain the missing optical flux. GRB 090417B 
is an example of a gamma-ray burst that is dark due to the localized dust structure in 
its host galaxy. 



Subject headings: dust, extinction - 
- gamma-ray burst: GRB090417B 



galaxies: individual (SDSS J135846.65+470104) 



Introduction 



The Swift observatory (JGehrels et al.ll2004h is a multi-instrument satellite mission that was 
designed to detect and rapidly localize gam ma-ray bursts (GRBs). The observatory contains three 
telescopes: the Burst Alert Telescope (BAT: iBarthelmv et al.l (|2005h ). the X-Ray Telescope (XRT; 
Burrows et all (J2005h . and the UltraViolet/Optical Telescope fUVOT: lRo~ming et al.l (12005^ . The 
BAT is used to identify GRBs and localize them to ~ 3' in the energy range 15-150 keV. Once the 
BAT has localized a burst Swift slews to point the XRT and the UVOT at the burst. The XRT 
obtains rapid X-ray localizations to < 5" in the energy range 0.2-10 keV while the UVOT obtains 
localizations to ~ / .'5 then repeatedly cycles through a set of optical and ultraviolet filters. Over 
the past five years Swift has detected GRBs at the rate of approximately 100 per year. Almost all 
of these Swift-detected GRBs have had X-ray afterglows, and optical or infrared afterglows have 
been detected for ~ 60% of them. The remaining « 40% of the Swift bursts have no reported 
optica l or infrared detections. In general between ~ 25-40% of GRBs are dark GRBs (jFvnbo et al 
2009bh . 



The gamma-ra y burst GRB 090417B was detected by the BAT as an image trigger at 13:17:23 
UT on 2009 Apr 17 (Isbarufatti et al.lbood v Swift immediately slewed to GRB 090417B and narrow- 
field observations began observing 387 s (XRT) and 378 s (UVOT) after the BAT trigger. The 
BAT light curve showed a long period of emission starting 200 s before the trigger (see Figure [1]). 

Several gr ound-based fo l low-up observations were made, but no optical or infrared afterglow 



was detected. iFvnbo et al.1 (|2009al ) found a source with R = 21.3 inside the XRT error circle 



(jSbarufatti et al.ll2009l ) and i dentified it as an SD SS galaxy (IStrauss et al.1 12002 ) and the possible 
host galaxy of GRB 090417B. iBerger Sz Foxl (|2009i ) measured a redshift of z = 0.345 for this galaxy 



based on the forbidden oxygen emission lines [O II] (3727) and [O III] (5006) as well as Ha emission. 
Dark GRBs have been an enduring mystery in the lore of GRB studies. In the pre-Swift era 



only about one qu arter of localized gamma-ray bursts had an optical afterglow (OA) detected (e.g., 



Fvnbo et al.ll200ll ). It was sometimes assumed that this was due to the response times of ground- 
based optical telescopes which allowed most afterglows to fade before they could be observed. 
Swift however, with its ability to train the XRT and UVOT on a burst within ~ 100 s of the 
initial detection, has been able to provide positions to within ~ 5" for most bursts that have been 
detected by the BAT. The rapid distribution of these positions has allowed ground-based optical 
and infrared follow-up observations for many of these bursts, often within minutes of the BAT 
trigger. Even with such rapid follow-up observations no optical or infrared afterglows have been 
found for ~ 40% of BAT-detected GRBs, thus demonstratin g that observational constraints were 



not the cause to the dark burst problem (JRoming et al.ll2006l ). 



Several attempts have been made to quantify the dark burst problem by establishing opera - 



tional definiti ons of dark bursts ( Jakobs son et al.1 12004 : iRol et al. 



2005 



van der Horst et al.1 12009 ) . 



Most recently Ivan der Horst et al.l ( 20091 ) proposed that dark GRBs be defined as those that have 
Pox < Px — 0.5 measured 11 hours after the burst. They define fiox as the spectral index between 
the optical and X-ray regimes and fix as the spectral index in the A-ray band. The spectral index 
is defined by f y {y) oc v~® where j v {y) is the flux density at a frequency v. This definition assumes 
that both the A-ray and optical afterglows are due to synchrotron radiation from a relativistic 
fireball. The synchrotron model predicts that the spectral slope will either be the same in both 
regimes (if there is no cooling break between them) or will differ by A/3 = 0.5 (if there is a cooling 
break between them). A difference of A/3 > 0.5 can not be explained in this model and thus 
requires that somethin g be suppressing the flux at optical wavelengths, resulting in a dark burst 
\ Jakobsson et al.l 12004 ) . 



Several possibilities have been proposed to explain dark bursts. First, several studies have 
suggested that dark bursts m ay be due to extinction along the line of sight, either in the Milky 



2008 



Way or in the host galaxy ( e.g., Piro et al.ll2002l ; lLevan et al.ll2006l : lBerger et al.ll2007l ; |jaunsen et al. 



Tanvir et al.ll2008l ). IPerlev et al.1 (J2009I ) suggest that this is the case for most dark GRBs. 



Further, A-ray observations show that some bursts with no op tical afterglows have systematically 
higher H I column densities than those with optical afterglows (IFynbo et al.ll2009bl ). This suggests 
that dark bursts may be suffering from higher extinction than optically-bright bursts. However, 
there is significant overlap in the distribution of Ah values between dark and bright bursts, so it is 
not clear if extinction is the sole parameter responsible for the darkness of these bursts. 



Schadv et al.l (J2007I . 120091 ) did panchromatic studies of several GRBs that were observed with 



Swift's XRT and UVOT (i.e., optically-bright GRBs) in order to probe the extinction in their 
host galaxies. They found that the combined A-ray /ultraviolet /opti cal data for most of these 
GRBs are best fit using a Small Magellanic Cloud (SMC) dust model (|Peilll992i ). The extinctions 
in these galaxies were small, with Ay < 1 mag. They also found a mean gas-to-dust ratio of 
Ah = (6.7 x 10 21 ) Ay, which is lower than, but s imilar to, th e value found in t he SMC (Ah = 
(15.4 x 10 21 ) Ay using equation (4) and Table 2 ofH dl992l ^. kann et al.l (hoOfih also found that 
the the gas-to-dust ratio for a large sample of host galaxies of optically-bright GRBs is lower than 



the SMC value. However, the gas-to-dust ratios in the host galaxies of dark GRBs tend to be 
more like that of the Milky Way than tha t of the SMC (e.g., 
Jaunsen et aPhoOsI : lElfasdottir et alJboogl ). 



Piroet al J 12002 : Berger et al.ll2007l : 



A second possible origin for dark bursts is that they are GRBs located at high redshift, so the 
Lyman break (or the start of the Lyman-alpha forest) is observ ed redward of the optical band. The 



most distant GRB observed to date is GRB 090423 at z = 8.2 (jSalvaterra et al 



2009 



Tanvir et 



ap 

2009), so some of the dark bursts could simply be at high redshift. However, iFvnbo et al.l (J2009b) 
have found that less than approximately 19% of GRBs have z > 7, whereas between 25% and 42% 
of GRBs are dark. Therefore, it is unlikely that all of the Swift dark burst lie at very high redshifts. 
The affects of relativistic beaming of GRB emission at very high redshift is poorly understood due 
to a lack of data on high-redsift bursts. We assume here that it does not affect our discussion in 
this paper. 

Third, the hypothesis that afterglow radiation is synchrotron radiation may be wrong. This 
could lead to a different spectrum from what is expected from the relativistic fireball model, in 
which case the above definition of a dark burst is invalid. 



Long-soft GRBs (JKouveliotou et al.lll993l ) tend to be found in host galax ies that are small 



i rregu lar, and have high specific star- formation rates (IChristensen et al.l 120041 ) . iLe Floc'h et al 
(|2003l ) found a median infrared luminosity of L ~ 0.08L* for the host galaxies of GRBs with optical 
afterglows and conclude that most hosts are sub-luminous at optical and near-infrared wavelengths. 
These GRB hosts tend to be small galaxies that have intrinsic luminosities similar to those of the 
dwarf galaxies in the Local Group. However, some hosts are significantly larger, such as that of 
GRB 990705 (JLe Floc'h et alJbooa ), which was an Sc spiral with L « 2L*. To date no long-soft 



GRB has been found in an elliptical galaxy. All of the hosts of long-soft GRBs have morpholo 



gies that are consistent with either exponential discs, or irregular structure (IConselice et al 



200 



5; 



Wainwright et al.ll2007l ). GRB hosts foll ow the size-luminosity r elation (e.g.. iTrujillo et al.ll2006l ) 
and extend it towards lower luminosities (IWainwright et al.ll2007l ). In spite of their generally small 
size the host galaxie s of long-soft GRBs tend to have high specific star- formation rates (SSFRs). 
Savaglio et al.l (2009j) found a median SSFR of 0.8 yr _1 , similar to what is seem in Lyman-break 
galaxies. 

It is important to realize, however, that these properties have been determined using GRBs 
that have usually been localized to sub-arcsecond precisions. In practice this usually means that 
the sample of GRB host galaxies is restricted to GR Bs that have had op tical afterglows. The host 
galaxies of dark GRBs may have different properties. IPerley et al.1 (J2009i ) identify host galaxies for 
14 dark GRBs and found that they have the same redshift distribution as the hosts of optically- 
bright GRBs. They also found that these hosts do not appear to be significantly different in size, 
structure, or luminosity from the hosts of GRBs with optical afterglows. They did find that a 
significant fraction of the hosts of dark bursts show evidence for high internal extinction [Ay > 2-5 
mag). However, many of their host associations are uncertain, with more than half of the galaxies 



having a probability of a chance association on the sky of P c h > 0.01. As we will show in this paper 
GRB 090417B is located 1'.'08 from the centre of a bright galaxy with a probability of a chance 
association of only P c h ~ 10 -3 . This makes GRB 090417B one of the strongest associations of a 
dark GRB with a host galaxy. 

In this paper we present space- and ground-based gamma-ray, X-ray, ultraviolet, and optical 
observations of GRB 09041 7B and its host galaxy. We have adopted a standard ACDM cosmology 
with a Hubble parameter of Hq = 71 km s' 1 Mpc -1 , a matter density of VL m = 0.27, and a 
cosmological constant of Q,\ = 1— £l m = 0.73. For this cosmology a redshift of z = 0.345 corresponds 
to a luminosity distance of 1.806 Gpc and a distance modulus of 41.29 mag. One arcsecond 
corresponds to 6.51 comoving kpc, or 4.84 proper kpc. The look back time is 3.80 Gyr. 



2. Data 

2.1. BAT Data 

The BAT scaled-map, event and DPH data were reduced using the standard BAT software 
available through HEASARC0. 

The T 90 duration of GRB 090417B was > 2130 ± 50 s in the 15-150 keV band. The light curve 
evolved gradually, and there were no features in the unimaged count rates that might indicate 
how long the burst persisted in the BAT energy range past that time. Even at that lower limit 



on T90, this was the longest GRB ever seen by BAT or BATSE (Fishman e t al.lll988r ). although 
GRB 060218 approached this duration with T 90 = 2100 ± 100 (|Oampana et allhoOfil ). 



The BAT light curve for GRB 090417B showed four broad, overlapping peaks in the flux at 
approximately T — 70, T + 400, T + 500, and T + 1600 s, as shown in Figure [TJ The 1-s peak flux of 
0.3 ±0.1 photons cm' 2 s _1 in the range 15-350 keV was approximately T + 500 s. The total fluence 
from T — 320 to T+2105 s was 8.20_ 21 x 10 -6 erg cm -2 in the 15-150 keV band. The spectrum was 
best fitted by a simple power-law function with an average photon index of 1.89 ± 0.12. The index 
did not change significantly during the burst. All the above uncertainties are at 90% confidence. 

The most impressive characteristic in the prompt emission of GRB 090417B was its extreme 
length. The burst was marginally above background (2.3a) on the 64-s interval starting immediately 
after a pre-planned slew maneuver ended at T — 320 seconds, although not at a high enough level 
to be detectable independent of the later emission. From that time, the flux increased steadily, to 
be detected as a point source on-board on the interval T — 64 to To s, and finally as a sufficiently 
significant source to trigger a burst response on the interval To to T + 320 s. Although it appeared 
to be declining toward background levels, the burst continued to be detected until the spacecraft 
slewed away due to the changing Earth constraint at T + 2105 s. 



http://swift.gsfc.nasa.gov/docs/software/lheasoft/download.html 
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Fig. 1. — The top four panels show the BAT light curves for the prompt emission from GRB 090417B 
in four energy bands. From the top down the energy bands are: 14-24 keV, 24-51.1 keV, 51.1— 
101.2 kev, and 101.2-194.9 keV. The bottom panel shows the total light curve in the 14-194.9 keV 
band. 



2.2. XRT Data 

2.2.1. Data Reduction 
The XRT on board Swift began o bservations on 2009 April 17 at 15:26:33 UT, 387 s after 



the BAT trigger (jSbarufatti et al.ll2009l ). and ended on 2009 May 5, with a total net exposure of 
1456 s in Windowed Timing mode (WT) and 83.36 ks in Photon Counting mode (PC). The XRT 
observations are listed in Table [TJ 

The XRT data were processed using the FTOOLS software package distributed inside HEA- 
SOFT (v6.6.3). We ran the task xrtpipeline applying calibrations and standard filtering criteria. 
Events with grades 0-2 and 0-12 were selected for WT and PC data respectively. The analysis was 
performed inside the 0.3-10 keV energy band. 

The best position for the afterglow was obtained using 6972 s of overlapping XRT Photon 
Counting mode data and 9 UVOT v-band images in order to correct the XRT astrometry mak- 
ing use_of_the XRT^UVOT_ahg2miejrt andniatching to the USNO-B1 catalogue, as described 



by iGoad et al.1 (|2007|) and Evans et al.1 (|2009h . Our best position is RA(J2000.0) = 209?6942 
(13:58:46.62), Dec(J2000.0) = +47?0182 (+47:01:05.4) with an error of l'/4 (90% confidence, in- 
cluding boresight uncertainties). 



2.2.2. Temporal Analysis 

Source photons were extracted from a region with a 30 pixel radius (1 pixel = 2'/36), with 
the exception of PC mode data for sequence 00349450000 where we used an annular region with 
radii 3 pixel and 30 pixel in order to correct for pile-up. The background was estimated from a 
circular region with a 50 pixel radius located away from any detected source in the field. When 
the count-rate dropped below the level of ~ 10~ 2 count s _1 we used the SOSTA tool of ximage, 
which corrects for vignetting, exposure variations and PSF losses within an optimized box, using 
the same background region. The 0.3-10 keV light curve thus obtained is shown in Figure [2] (top 
panel). The light curve has been rebinned in order to achieve a minimum signal-to- noise ratio 
(S/N) of 3 for each point. The best fit to the light curve is given by a doubly broken power law 
with indices a x = 0.86 ± 0.03, a 2 = 1.40 ± 0.04, and a 3 = 2.0tol, and breaks at T + 9100+ioq s 
and T + 4.5_ 1 ' 2 x 10 5 s. The second break could be interpreted as a jet break. However the lack of 
an optical detection for this burst does not allow us to confirm or deny this hypothesis. Also see 
§ 14.31 for an alternate interpretation of the late-time steepening of the X-ray light curve. 

The first sequence of X-ray data is dominated by a strong flare. There is a rise that began 
before the start of XRT observations and peaked around T + 530 s. This is followed by a double- 
peaked flare which started at T + 1260 s, peaked at T + 1410 s and T + 1470 s, and had a decay 
time of 240 s. Further flaring activity is apparent up to at least T + 4 x 10 5 s. The bottom panel 
of Figure [2] shows the hardness ratio for the afterglow measured using the 0.3-2 and the 2-10 keV 



bands. The first orbit (up to T + 1300) shows a hard hardness ratio with spectral evolution that 
follows the flaring activity. The same is true for the remaining observations, with the hardness 
ratio slowly decreasing with time. 



2.2.3. Spectral Analysis 

In order to detect possible spectral variations we extracted spectra separately for the two flares, 
the early decay part from T+5000 s to T+9100 s, the steep decay from T+9100 s to T+4.5 x 10 5 s, 
and the late decay from T + 4.5 x 10 5 s to the end of observations. Source and background spectra 
were extracted using the same regions used for the light curves, with the exception of the late 
decay spectrum where the count rate was too low to group the energy channels in order to have 
a Gaussian distribution of the number of photons per grouped channel. Hence, for this spectrum 



we replaced the x 2 statistic with the C statistic (|Cashlll979l ). which can be used whenever the 
Gaussian approximation is not valid, provided that the contamination from background photons is 
negligible. To ensure this, we extracted photons from a region with a 10 pixel radius. 

The ancillary response files were generated using the task xrtmkarf, and with the excep- 
tion of the late decay spectrum we grouped channels in order to have at least 20 photons per 
bin. Spectral fitting was performed using XSpec (vll.3.2). For each spectrum we modelled the 
Galactic absorption u sing a warm absorber with neutral hydrogen column density 1.60 x 10 21 cm -2 



(JKalberla et al.ll2005l ) and included an intrinsic warm absorber at the redshift of the host galaxy 
(z = 0.345). The best model for the spectra of the flares was a power law with a high energy cutoff, 
while for the remaining parts of the afterglow the best fit was obtained using a single power law. 
The best fit parameters are given in Table [2j 

The fact that the first phases of the XRT observations are best fitted by a cutoff power law 
with the cutoff energy around 3 keV indicates that the peak energy of the flares was passing through 
the XRT band at the time. After the flares the spectrum is described by a power law with photon 
index of 2.0 ±0.1. After the first break, up to T + 29.6 ks, the spectrum is described by a power law 
with similar parameters (r = 2.3 ± 0.2 and an absorbing column Nr = 2.0 ± 0.4 22 cm -2 ). At later 
times the spectrum softens significantly to a photon index of T = 3.8 ±0.5 and shows no significant 
variations across the break at T + 4.5 x 10 5 s up to the end of the observations. 

The intrinsic absorbing column is in the range (1.1 — 2.4) x 10 22 cm -2 . From our analysis of 
the spectral energy distribution (§ [372]) we find that the best fit to the data is given by a Milky Way 
extinction law instead of the SMC-like extinction law found in many GR B host galaxies. Therefore , 



we use the relationship between JVh and extinction in the Milky Way from lPredehl fc Schmitti (119951 ) 
to find that the observed absorbing column corresponds to an extinction A v « 11 mag. This high 
value for the extinction in the host along the line of sight to the burst is capable of explaining the 
dark optical nature of this event. 



Table 1. Journal of XRT observations. 



Sequence 


Start Time 


End Time 


Start time 1 


End Time 1 


Obs. Mode 2 


Exposure 




UT 


UT 


(s) 




(?) 






(<0 


00349450000 


2009-04-17 15:26:40 


2009-04-17 23:29:46 


3.932 


X 


10 2 


2.055 


X 


10 3 


WT 


1.456 x 10 3 


00349450000 


2009-04-17 15:38:43 


2009-04-17 23:33:26 


1.116 


X 


10 3 


2.960 


X 


10 4 


PC 


9.214 x 10 3 


00349450002 


2009-04-18 10:29:48 


2009-04-18 11:07:52 


6.898 


X 


10 4 


7.124 


X 


10 4 


PC 


2.252 x 10 3 


00349450003 


2009-04-18 16:57:16 


2009-04-18 20:20:44 


9.223 


X 


10 4 


1.044 


X 


10 5 


PC 


4.871 x 10 3 


00349450006 


2009-04-19 01:17:02 


2009-04-19 12:50:57 


1.222 


X 


10° 


1.638 


X 


10° 


PC 


9.709 x 10 3 


00349450007 


2009-04-20 04:41:58 


2009-04-20 22:26:57 


2.209 


X 


10° 


2.848 


X 


10° 


PC 


6.131 x 10 3 


00349450008 


2009-04-21 14:13:55 


2009-04-21 23:59:56 


3.416 


X 


10° 


3.768 


X 


10° 


PC 


7.910 x 10 3 


00349450009 


2009-04-22 04:29:42 


2009-04-22 14:20:47 


3.930 


X 


10° 


4.284 


X 


10° 


PC 


4.574 x 10 3 


00349450010 


2009-04-23 00:07:52 


2009-04-23 22:40:57 


4.637 


X 


10° 


5.448 


X 


10 5 


PC 


2.989 x 10 3 


00349450011 


2009-04-24 00:18:14 


2009-04-24 22:57:56 


5.507 


X 


10° 


6.323 


X 


10° 


PC 


1.344 x 10 3 


00349450012 


2009-04-25 11:39:16 


2009-04-25 23:01:57 


6.779 


X 


10° 


7.189 


X 


10° 


PC 


1.480 x 10 3 


00349450013 


2009-04-26 02:11:18 


2009-04-27 23:12:56 


7.303 


X 


10° 


8.923 


X 


10 : ' 


PC 


2.140 x 10 3 


00349450014 


2009-04-27 00:40:17 


2009-04-27 21:39:56 


8.112 


X 


10° 


8.868 


X 


10° 


PC 


1.224 x 10 3 


00349450015 


2009-04-28 00:37:57 


2009-04-28 23:20:56 


8.975 


X 


10° 


9.792 


X 


10° 


PC 


8.368 x 10 3 


00349450016 


2009-04-29 00:41:42 


2009-04-29 23:26:57 


9.841 


X 


10° 


1.066 


X 


10° 


PC 


4.537 x 10 3 


00349450017 


2009-05-01 00:40:42 


2009-05-02 23:37:57 


1.157 


X 


10° 


1.326 


X 


10° 


PC 


1.604 x 10 4 



'Time since BAT trigger. 

2 WT, window timing mode; PC, photon counting mode 



Table 2. Spectral parameters. 



Start time 
(<0 



End Time 1 
(<0 



Mi 



(10 2 



cm 



Cutoff Energy 
(kcV) 



Mean Flux 
(erg cm" 2 s~ 4 ) 



X 2 /d.o.f. 



3.932 x 10 2 
1.116 x 10 3 
5.000 x 10 3 
9.100 x 10 3 
6.898 x 10 4 
4.637 x 10 5 



1.116 x 10 3 
2.055 x 10 3 
9.100 x 10 3 
2.960 x 10 
4.284 x 10 5 
1.326 x 10 6 



1.5 ±0.2 
1.3 ±0.2 
2.0 ±0.4 

1.6 ±0.4 
1.9 ±0.5 
< 1.9 



0.2 ±0.2 
0.3 ±0.2 
1.9 ±0.2 
2.3 ±0.2 
3.8 ±0.5 
3.1 ±0.9 



3.0 

2.9 



+0.7 
-0.5 
+0.5 

0.4 



1.3 x 10~ 09 
2.8 x 10~ 09 
6.2 x 10~ n 
1.7 x 10" 11 

4.4 x 10" 13 
7.7 x 10~ 14 



446/421 

604/523 

39/45 

15/29 

13/22 

32/24 2 



1 Time since BAT trigger. 
2 C-statistic/bins 
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Fig. 2. — Top: The X-ray flux light curve of GRB 090417B. Bottom: The hardness ratio evolution 
of GRB 090417B in the 0.3-2 and 2-10 keV bands. The blue points are window timing mode data 
and the red points are photon counting mode data. 



11 



2.3. UVOT Data 



The Swift/UYOT began observations 378 s after the BAT trigger (ISbarufatti et al.ll2009l ) with 
a 9 s settling mode exposure using the v filter. This was followed by a sequence that rotated through 
all seven of UVOT's lenticular filters. The first of these exposures was a 147 s exposure starting at 
395 s taken with the white filter. Observations continued until 1 325 874 s (= 15.3 days) after the 
trigger using primarily the u and ultraviolet filters. No optical or ultraviolet afterglow was detected 
in any of the UVOT data. Once it was determined that there was no optical or ultraviolet afterglow 
for GRB 090417B UVOT observed primarily with its ultraviolet filters in order to obtain ultraviolet 
magnitudes for the nearby galaxy that can not be obtained from ground-based observatories. 



The SDSS galaxy noted bv lFvnbo et al.l (J2009al ) is detected when we coadd the UVOT images. 
We determined a centroid for the SDSS galaxy of RA, Dec = 13:58:46.66, +47:01:04.4 (J2000.0) with 
an estimated internal uncertai nty of $'.§7 and an es timated systematic uncertainty relative to the 
USNO-B1.0 catalogue of 0?42 (jBreeveld et al.ll201ol ). These uncertainties are the 90% confidence 
intervals. This galaxy is l'/08 southeast of the centre of the UVOT-enhanced XRT error circle. The 
field of GRB 090417B is shown in Figure El The SDSS galaxy is well-isolated from other sources in 
the field, so there is no contamination from neighbouring sources when doing aperture photometry. 

We obtained the UVOT data from the Swift Data Archive on 12 May 2009. This data has 
had bad pixels identified, mod-8 noise corrected, and has been transformed into FK5 coordinates. 
We used the standard UVOT data analysis software distributed with HEASOFT 6.6.2 along with 
the standard calibration data. Photometry was done on the SDSS galaxy using UVOTSOURCE with 
a circular source aperture of radius 4'.'0 for the source region and a circular aperture of radius 15" 
centred at RA, Dec = 13:58:47.57, +47:00:33.3 (J2000.0) for a background region. The background 
region was selected to have similar background properties to those at the location of the galaxy, 
and to be free of contaminating sources. The UVOT photometry of the SDSS galaxy is presented 
in Table [3l We find no evidence for any change in the magnitude of the galaxy, in any filter, during 
the course of the U VOT observations. Our photometry was calibrated to the UVOT photometric 
system described in iPoole et al.l (|2008l ) . 



The UVOT-enhanced XRT position for GRB 090417B is RA = 13:58:46.62, Dec = +47:01:05.4 
(J2000.0) which corresponds to Galactic coordinates of £ u ,b u = 93? 7 486, +66? 1127. The line-of- 
sight Galactic reddening in this direction is Eb~v = 0.02 ± 0.01 mag (jSchlegel et al.lll998l ). This 
corresponds to extinctions in the UVOT filters of of A v = 0.05, Aj, = 0.07, A u = 0.08, A uvw \ = 0.11, 
Auvm2 = 0.16, and A uvw2 = 0.14, and A white = 0.08 mag. 

Since there is a candidate host galaxy for GRB 090417B it is not possible to put a direct faint 
limit on th e luminosity o f the o ptical afterglow beyond the observed magnitude of the galaxy itself 
(R = 21.3 iFynbo et al.l (j2009al )). In order to constrain the maximum possible luminosity of the 
afterglow we used the UVOT data to look for variability in the galaxy's magnitude. To do this we 
divided the data into early and late epochs. The early epoch consisted of all the UVOT data taken 
during Swift's initial automated observations of GRB 090417B (OBSID 00349450000), which covers 
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Fig. 3. — This Figure shows the averaged NOT/ALFOSC i?-band image of the field of 
GRB 090417B with the XRT error circle. The SDSS galaxy is in the southeast quadrant of the 
XRT error circle. The arms of the compass have a length of ten arcseconds. 

Table 3. This Table lists the UVOT photometry of SDSS J135846.65+470104.5, the host galaxy 
of GRB 0904 17B. Exposure is the total exposure time in each filter. All upper limits are 3-<r 

upper limits. 

Filter Exposure Mag Err 



V 


1237 


>20.9 




b 


862 


>21.4 




u 


22 554 


23.09 


0.38 


uvwl 


10 850 


23.12 


0.56 


uvm2 


15 942 


22.49 


0.29 


uvw2 


30 827 


22.69 


0.19 


white 


841 


>22.2 
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times up to 29 604 s (~ 8 hr) after the BAT trigger. This is the period when the optical afterglow 
is expected to be brightest. The late epoch included all observations later than this. The UVOT 
exposures for each epoch were coadded to produce two deep images for each UVOT filter. We 
photometered the galaxy in each image and recorded the difference between the two magnitudes. 
If we assume that the first epoch magnitude is due to the host+afterglow and the second epoch 
magnitude is due to the host alone then the difference between these two magnitudes puts an upper 
limits on the luminosity of the afterglow. We find no evidence for any change in the magnitude of 
the host galaxy between the first and last epochs larger than ±0.2 mag. Therefore, the afterglow 
could not have contributed more light than what corresponds to a 0.2 mag increase in luminsity. 
From this we obtain the following constraints on the magnitude of the optical afterglow at 4 hr (= 
14 800 s) after the BAT trigger: u oa > 24.9, uvwl oa > 24.9, uvm2 oa > 24.1, and uvw2 oa > 24.4. 
The uncertainty in these estimates is about 1.0 mag. 



2.4. Nordic Optical Telescope Data 

We obtained optical images of GRB 090417B on 2009 Apr 17-18 with the ALFOSC-FASU 
mounted on the Nordic Optical Telescope (NOT). ALFOSC-FASU was operated in its high gain 
state (0.726 e _ /ADU). The read-out noise was 3.2 e _ /pixel and the pixel scale was O'.'19/pixel. The 
NOT i?-band image is shown in FigureEJ Data reduction was done following the standard procedure 
for optical CCD data and the individual images were aligned by centroiding on several stars in 
each exposure. Photom etry was performed on the combined NOT images using SExtractor 
( Bertin &: Arnoutslll996l ). We calibrated the data using the standard ALFOSC zero pointqj. We 



find B = 23.25 ± 0.13 (Bessel B),R = 21.34 ± 0.03 (Bessel R), and % = 20.82 ± 0.07 (interference 
i) for the SPSS gala xy. The Galactic extinctions for the NOT data are Ab = 0.07, An = 0.05, 
(jSchlegel et allllflfld ). and A t = O.ofl 



3. The Host Galaxy 
3.1. Probability that the SDSS Galaxy is the Host 



We used the methodology of iBloom et al.1 (|2002l ) to estimate the probability that a galaxy 



would lie within l'/08 of the centre of the XRT error circle by chance, P&. This assumes that 
galaxies are randomly distributed on the sky and that there is no clustering. This assumption is 
not formally correct as galaxies do cluster. However, an examination of the field of GRB 090417B 
suggests that this GRB did not occur in a significant over-density of background galaxies, so 



2 http://www.not. iac.es/instruments/alfosc/zpmon/ 
3 http://irsa.ipac. caltech.edu/applications/DUST/ 
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distortions in the computed P r ^ va lue are likely to be small. iBloom et al.l (J2002I ) use the galaxy- 
number counts of lHogg et al.l (|1997l ) to generate the surface density of galaxies with R magnitudes 
brighter than some limiting magnitude. In this case the R band is for the Keck 10-m Low Resolution 
Imaging Spectrometer (L RIS) R filter. We a dopt R = 21.34 ± 0.03 (§ I2.4p as the magnitude of the 
SDSS galaxy. Using the IBloom et al.l ((2.002) formula for an error circle containing a galaxy this 
yields P ch = 0.003. This is a fairly small value, which suggests that SDSS J135846.65+470104.5 is 
likely to be the host galaxy of GRB 090417B. 

A chance alignment probability of P c h ~ 10~ 3 is large enough that one can not rule out the 
possibility of a misidentification. A visual examination of the NOT images shows that the XRT error 
circle includes a significant fraction of the SDSS galaxy (see Figure [3]). The observed separation 
between the X-ray afterglow and the centre of the galaxy is l'/08. At a redshift of z = 0.345 this 
corresponds to a projected separation of 5.23 proper kpc. The galaxy appears point-like in all of 
our data. The NOT i?-band data has a seeing FWHM of / /91, which repr esents an upper limit on 
the effective radius of the galaxy of r e < 4.4 proper kpc lBerger et al.1 (120091 ) found g' = 22.47±0.16 
for the galaxy, which implies an absolute magnitude of M y ~ —18.8. Thi s puts the galaxy inside 
the observed range of GRB host sizes and luminosities (see iFruchter et al.l (|2006l ) , their Fig. 4) . 



In order to test our hypothesis that SDSS J135846.65+470104.5 is the host galaxy of GRB 090417B 
we obtained the Chandra observations of the X-ray afterglow that were taken on 2009 May 11, 24 
days (~ 2 x 10 6 s) after the burst (Proposal Number: 10900117, PI: Andrew Levan). An X-ray 
source (consistent with the XRT position) is clearly detected at the position of the SDSS galaxy. 
The observed flux from this source is consistent with the X-ray flux that is predicted from the XRT 
light curve assuming a power-law decay with 03 = 2. This suggests that the Chandra source is 
the X-ray afterglow of GRB 090417B, and that there is minimal (if any) contamination from the 
galaxy. 

Taken together these arguments strengthen the case that SDSS J135846. 65+470104. 5 is the 
host galaxy. Therefore, we conclude that we have identified the host of the dark burst GRB 090417B 
to a high degree of confidence. 



3.2. Spectral Energy Distribution 



We used the combined ultraviolet, optical, and infrared photometry of the SDSS galaxy to 
determine its spectral energy distribution (SED). Photometry for this galaxy was taken from our 
observations, the SDSS Web sitqj, and the various GCN Circulars on this burst, and is listed 
in Table HI We used XSpec to fit a power law to the observed SED. The best fit, with no 
host extinction, has a power-law index of T = 2.67_ 'o3- However, the goodness of fit is only 
X 2 /dof = 81.40/12, which is very poor. Therefore, we also tried including extinction due to dust in 



1 http://www.sdss.org/DR7/ 
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the host galaxy when fitting the SED. Dust laws for the Small Magellanic Cloud, Large Magellanic 
Cloud and the Milky Way were tried. These led to a significant improvement in the goodness of 
fit with the best fit being found for Milky Way extinction in the host galaxy. Our best fit has a 
spectral index of T ~ 0.80^ 



~ ' 55 , and a host extinction of Ay 



o r + 1.0 
A -°-0.5- 



3.3. Luminosity and Star-Formation Rate 

SDSS J135846. 65+470104.5 has a colour of B-R = +1.91 ±0.13, making it a fairly red galaxy. 
The typical absolute magnitude for a red galaxy at 0.2 < z < 0.5 is (M B ) AB = —20.44 (JLillv et al 
19951 ) . assuming a cosmology with (i?o,^m>^A) = (50,1,0). For our adopted cosmology that 
corresponds to (M B ) AB = —20. At z = 0.345 the rest-frame B band corresponds roughly to the 
observed r' band. Assuming that the galaxy has a power-law spectrum (see § I3.2p and r' = 21.6, 



then it has a rest-frame luminosity in the B band of Lb ~ 1.3L^ where L B is the rest-frame -B-band 
luminosity of a typical red galaxy at z = 0.345. This indicates that the host is approximately an 
L* galaxy with a total luminosity that is roughly similar to that of the Milky Way. However, we 
stress that this result is somewhat uncertain because M* is highly corre lated with the slo pe of the 
faint end of the galaxy luminosity function, and with its normalization (jLilly et al.lll995l ). 



We can estimate the star- formation rate in the host galaxy using Eq. 2 of lMadau et al.l (jl998l ). 
At the redshift of the burst 1500 A approximately corresponds to the UVOT uvw2 filter while 
2800 A approximately corresponds to the SDSS v! filter. Correcting for Galactic extinction these 
magnitudes for the host become uvw2o = 22.55 and u' r 

(in M Q yr" 1 ) of 0.3 (1500 A) and 0.7 (2800 A) assmnuiR a. 



n - 2 3.15. These yie ld star- formation rates 



Salpeterl (| 19551 ) initial mass function, 



and 0.6 (1500 A) or 1.2 (2800 A) assuming a IScalol ( 19861 ) initial mass function. These values 
assume that there is no extinction in the host galaxy. This is clearly not correct, so they actually 
represent a lower limit on the star- formation rate of approximately 1 A4q yr _1 . Correcting for 
the mean extinction in the host (A v = 3.5 mag) gives a star-formation rate of > 10 2 M@ yr _1 . 
However, this corrected values should be treated with caution since extinction in the host is likely 
to be highly variable and may not be affecting all star-forming regions equally. 

The host galaxy appears to be roughly similar in luminosity to the Milky Way, but with a 
significantly higher star-formation rate. This puts it on the massive end of the GRB host galaxy 
distribution. 



4. The Gamma-Ray Burst and its Afterglow 
4.1. Is GRB 090417B a Dark Burst? 



No afterglow has been detected for GRB 090417B at infrared, optical, or ultraviolet wave- 
lengths. Changes in the near-ultraviolet flux from the host galaxy between early and late times are 
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Table 4. This Table lists all of the available photometry of SDSS J135846.65+470104.5, the host 

galaxy of GRB 090417B. 

Filter Mag Err Source 



uvw2 


22.69 


0.19 


1 


uvm2 


22.49 


0.29 


1 


uvwl 


23.12 


0.56 


1 


"UTJVOT 


23.09 


0.38 


1 


v! 


23.23 


0.53 


2 


^UVOT 


>21.4 




1 


B 


23.25 


0.13 


1 


g' 


22.02 


0.16 


2 


g' 


22.47 


0.16 


3 


"WTJVOT 


>20.9 




1 


R 


21.34 


0.03 


1 


r' 


21.9 


0.3 


4 


r' 


21.62 


0.09 


2 


r> 


21.62 


0.10 


3 


i 


20.82 


0.07 


1 


i' 


21.41 


0.11 


2 


i' 


21.31 


0.12 


3 


z' 


20.78 


0.25 


2 


z' 


21.49 


0.30 


3 


J 


20.0 




5 


K s 


18.5 




5 
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consistent with no contribution from an afterglow down to a limiting magnitude of u oa > 24.9. We 
are assuming that the bright limit on the afterglow magnitude is valid for all time after four hours 
after the BAT trigger. There is no evidence for a brightening of the host galaxy at late times, so 
this is a reasonable assumption. Therefore we adopt this as an upper limit to the luminosity of the 
afterglow at 1 1 hr and compare it t o the X-ray flux at 11 hr. This yields fiox < —1.9, making it 
a dark burst (jJakobsson et al.ll2004l ). T he X-ray spectrum at this time is evolving from fix = 1.3 
to fix = 2.3, so using the definition of Ivan der Horst et al.l (J2009J) we find that fiox < fix — 0.5, 
also making GRB 090417B a dark burst. Our limit on fiox is somewhat uncertain due to the large 
uncertainty in ou r estimate of tt m . How ever, for GRB 090417B to be an optically bright burst by 
the definition of (JJakobsson et al.l 120041 ) then the optical afterglow would have had to have had 
u < 21.7, which is inconsistent with the observed host+afterglow magnitude (u = 23.09 ± 0.38). 



Our computation of the spectral index between the optical and X-ray regimes depends on the 
limiting magnitude of the afterglow in the UVOT u band as determined from the statistical error 
in the observed magnitude of the host galaxy at two epochs. In order to test that GRB 090417B 
really is a dark burst we compute the expected UVOT u-band magnitude for the afterglow given 
the observed X-ray spectrum and assuming that the cooling frequency is between X-ray and optical 
wavelengths. This yields a prediction of u < 21.7 at 11 hours, which is ruled out by the observations. 
Therefore, we conclude that GRB 090417B is a dark GRB. 

The difference between the predicted u magnitude and the upper limit that we derive in § 12.31 
is Au > 3.2 mag. If we assume a Milky Way extinction law this corresponds to Ay > 2.5 mag of 
extinction in the host galaxy along the line of sight to the burst. 



4.2. Closure Relations 



If we assume the canonical X-ray light curve for GRB afterglows ( Nousek et al.l l2006; Z hang et al 



20061 ) then we can use the temporal decay indices and the spectral slopes to constrain the underly- 
ing physics. We stress that the following analysis is only valid if we assume the relativistic fireball 
model. Alt ernate models for the X -r ay afterglows of G RBs are available, such as the cannon- 
ball m odel (|Dar fc De RuiulalhoOOal lbl: bado et al.ll2002alf bh. th e disc model of iGannizzo fc Gehrels 



model (JYamazakil 120091 ; iLiang et al. 



i II -| |— " 

reverse shock model (IGenet et al,ll2007l: lUhm fc Beloborodovll2007l ). and the prior outflow 



200 1 



9) 



( 2009J ). the late inter nal emission mode l (Ghisell ini et al. 2007 ; Kumar et al.ll2008l ). the long-lived 



emission 



The time and post-break decay index of the late-time break in the X-ray light curve at 



4.5_ 1 ' 2 x 10 5 s after the BAT trigger are consistent with a jet break. However, the dust scat- 
tering interpretation of the late-time spectral softening (see § 14. 3p suggests that the steepening of 
the light curve after ~ 4.5 x 10 5 s is due to scattered emission due to dust dominating over the syn- 
chrotron component of the afterglow. Unfortunately there is no data in other wavelength regimes 
to test if the late-time steepening is achromatic (which would argue for a jet break), therefore it is 
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uncertain if this is a jet break. 



Sari et al.l (|1999l ) give the relationships between the temporal and spectral slopes and the 



elec tron power-law distrib ution index, p, for a blast expanding into a homogeneous external medium, 



and I Chevalier &: Lil ( 19991 ) give the relationships for a stellar wind (p(r) oc r 2 ) envir onment. Both 



of the se papers assume that p > 2. For th e case were p < 2 the relationships of iDai fc Cheng 
(|200ll ) are used. iPanaitescu &: Kumarl (|2002l ) found that GRB afterglows exhibit a range of electron 
distributions with values for the ten GRBs they studied being between p = 1.4 and p = 2.8. 

The X-ray spectral index during the first part of phase III of the afterglow (9100 < t < 29 600 s) 
is fix = T — 1 = 1.3 ± 0.2. If the cooling frequency is above the X-ray regime during this time 
then the predicted electron index i s p = 20x + 1 = 2(1 . 3 ± .2) + 1 = 3.6 ± 0.4, which is well 



outside the range of values found by IPanaitescu &: Kumarl (120021 ) . If the cooling frequency is below 
the X-r ay band during this time the n p = 2(3x = 2(1.3 ± 0.2) = 2.6 ± 0.4, which is consistent 
with the IPanaitescu &: Kumarl (J2002I ) range of electron indices. For v c < i>x the spectral index 
predicts a temporal decay of a x = 3/2/3 x - 1/2 = 3/2(1.3 ± 0.2) - 1/2 = 1.45 ± 0.30 for either a 
homogeneous interstellar medium (ISM) or a stellar wind environment. The observed X-ray decay 
during this time is ax = 1-40 ± 0.04, which is consistent with the predicted value. Therefore, we 
conclude that the cooling frequency is most likely below the X-ray band at 9100 < t < 29 600 s. 
Unfortunately it is not possible to tell, from the X-ray data, if the burst is expanding into a uniform 
external medium or a pre-existing stellar wind. The broad agreement between oli and fa with the 
predictions of the fireball model suggests that this is a reasonable explanation for the physics of 
this afterglow during this phase of the afterglow. 

Between 68 980 s and 428 400 s the spectral index softens to 2.8 ± 0.5. In this interval the 
fireball model predicts ax = 3.7 ± 1.0 for a homogeneous ISM and a wind if u c < vx- If v c > ^x 
then the environment must be windy (since v c is increasing) and ax is predicted to be 4.7 ± 1.0. 
The observed value is 02 = 1.40±0.04. Neither case is consistent with the data in this time interval. 

The data before 68 980 s is broadly consistent with a relativistic fireball expanding into either 
a homogeneous ISM or a wind-stratified external medium, and a cooling frequency below the X-ray 
band in the period 9 100 < t < 29 600 s. However, this model can not explain the data later than 
68 980 s. Some other physical process is needed to explain the late-time softening of the X-ray 
spectrum. 



4.3. Dust Scattering 

For GRBs occurring in dusty star forming regions, the X-rays emitted during the prompt 
phase may be scattered by dust grains i n small angles, which leads to a delayed X-ray afterglow 
component (JMeszaros &; Gruzinovl l2000l ) . Placing circumburst dust at ~ 10-100 pc away from 
the burst, the calculated dust echo X-ray lightcurve displays a shallow decay l ightcurve followe d 
by late steepening, as is commonly observed in many GRBs (jShao fc Daill2007l : IShao et al.ll2008l ). 
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One major featur e of this model is the strong spectral softening evolution in the echo emission 
(jShen et al.ll2009l ). Such a softening feature has not commonly been observed in the majority 
of GRBs, so the dust scattering model is not favored to interpret those bursts. However, GRB 
090417B shows a spectral softening during the shallow decay phase. Also the optically dark nature 
of the afterglow suggests that there is a substantial amount of dust in the circumburst environment. 
All these make this burst a very probable candidate for the dust scattering model. In modeling this 
burst, we assume a dust sh eet with dis t ance Rd from the burst, with other dust grain properties 
the same as in Figure 1 of IShen et al.l (|2009i ). We find the calculated light curve can reproduce 
quite well the XRT light curve from T — Tq = 5 x 10 3 s to 10 5 s (we consider the flares prior to this 
time as not being due to dust scattering because the 7-ray emission was still on at that time) for 
R d rs 30-80 pc and for a broad range of the source emission spectral index values. 

A strong test of the model is whether the spectral index evolution can be reproduced. In 
Figure H] we plot the temporal evolution of the observed XR T spectral index. T he dust scattered 
emission spectral index is calculated in the same way as in IShen et al.l ()2009|). The amount of 
softening predicted by the model is large, i.e., Afi ss 3-4, and it is independent of the source 
(prompt) emission spectral index. We tried to use the reported prompt BAT spectral index, fin AT = 
0.89, and the even harder spectral indices observed in the two X-ray flares (T — Tq = 400 — 2 x 10 3 s), 
fix = — 1) as the source X-ray spectral index, but found that the model predicted emission is 
significantly softer than what is observed. However, if one assumes an even harder spectral index of 
fix = —2 for the prompt X-rays, both the X-ray light curve and spectral evolution can be properly 
reproduced (see Figured]). Such a hard X-ray spectrum is consistent with the self-absorbed regime 
of a synchrotron or inverse Compton scattering spectrum. We therefore assume that the GRB 
prompt emission spectrum has a break above the X-ra y band due to the sel f-absorption frequency, 
u a , below which the spectrum has fi — —2 or —5/2 fe. g.JGranot &: Sarill2002l ). A large v a is possible 

The best model fit to 



Shen&ZhanE 



Jl2Q09|). 



if the prompt emission radius is small enough (e.g. 

the lightcurve and the spectral evolution is shown in Figured! The fit to the X-ray light curve has 

a reduced \ 2 of 2.7, and the fit to the spectral index evolution has a reduced \ 2 of 1.5. 

We also estimate the amount of dust in terms of rn , the scattering optical depth to the 1 keV 
photons, using the same formula as in I S hen et al.l (2009j). The source X-ray fluence is extrapolated 
from the observed 7-ray fluence using the assumed spectral form, and the scattered emission X-ray 
fluence is estimated from the observed XRT fluence for T — Tq > 5 x 10 3 s using an intermediate 
value of observed spectral index fix = 1.5. We find To ~ 2.5. This means that over half of prompt 
1 keV photons would be scattered by the dust. This tq value also implies A y sa 15 — 40 mag , 



according to some empirical relations between the two dust properties used in IShen et al.l ((2009). 
This Ay is much higher than the average value for the host galaxy that was derived from the 
spectral energy distribution in Section 13.21 but is consistent with the extinction derived from the 
X-ray measurement of the hydrogen column along the line of sight to the burst, Ay ~ 11 mag. 
This high line-of-sight extinction naturally explains the non-detection of an optical afterglow for 
GRB 090417B. Our model indicates that the dust grain column density along the line of sight is 
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« 10 12 cm 2 , and the density of the intervening dust sheet is n > 10 8 cm 3 . 

The data of GRB 090417B is consistent with the dust scattering model under the condition 
that the prompt X-rays have a self-absorbed spectrum in the X-ray range. The large amount of 
dust inferred from the X-ray modeling is consistent with the dark nature of the optical afterglow. 



4.4. Energetics of the Burst 

BAT observations of GRB 090417B were cut off at approximately 2105 s after the BAT trigger 
when Swift went into an Earth limb constraint. Therefore the total duration, and thus the total 
fluence, of the burst is not known. However, the fluence while the burst was observed was 8.20l 21 x 
10 -6 erg cm -2 in the 15-150 keV band. This puts a lower limit on the total fluence. We can use 
this lower limit to estimate a lower limit on the isotropic equivalent energy of GRB 090417B. The 
BAT spectrum is best fit by a single power law, so at a redshift of z = 0.345 this corresponds to an 
isotropic equiv alent energy of E\ Kn > 2.4_ ' 6 x 10 51 erg. We performed a k correction following the 



prescription of iBloom et al.1 (|200ll ) and find k = 2.6 (with some uncertainty due to the arbitrary 
choice of cut-off energy for the power-law spectrum). This yields a ^-corrected total isotropic energy 
of -Ei S o > 6.3 x 10 51 erg in the 20-2000 keV band. Our data does not allow us to constrain the 
time of a jet break, so we are unable to estimate the jet opening angle or the intrinsic gamma-ray 
energy of GRB 090417B. 



4.5. Constraints on a Supernova Component 

To date no supernova signature has been seen for any dark GRB. If a supernova were detected 
in a dark GRB that would have profound implications for origin of that dark burst. The detection of 
a supernova would indicate that the burst was not dark due to extinction or high redshift, but must 
have a fundamentally different afterglow than conventional optically-bright bursts. In light of this 
we examined our data to see if there is any evidence for a supernova component for GRB 090417B. 

UVOT observations were taken up to 15 days after the BAT trigger, so we searched this 
late-time data for evide nce of a supernova component. A type Ib/c supernova like SN1998bw 



(jPatat &: Piemontel ll998) at a redshift of z = 0.345 is expected to peak at approximately 10(1+2) ~ 
13 days after the burst. We find no evidence for a change in the luminosity of the host galaxy down 
to Au ~ 0.2 mag, which corresponds to a upper limit on the magnitude of the supernova of 
un m > 24.9. At z = 0.345 the central wavelength of the UVOT u band corresponds to Ao ~ 2600 A. 
No observational data is available for the peak magnitude of S N1998bw at 2600 A. However 



Swift /UVOT ultraviolet light curves are available for SN2007Y (jBrown et al||2009| ), a SN Ib/c 
with a peak absolute magnitude of uvwl ~ —17 in the rest frame. This corresponds to a rest frame 
uvwl of ~ 24 at the redshift of GRB 090417B. This is similar to our limiting magnitude, so if we 
assume that SN2007Y is typical of the SNe Ib/c associated with GRBs we are not able to constrain 
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Fig. 4. — Dust scattering model fit to GRB 090417B X-ray data. Top panel: the fit to the light 
curve. Bottom panel: the fit to the spectral index evolution. The model predictions are shown as 
solid lines and data as crosses. The model parameter values are same for both panels, and are as 
labeled. The source spectrum is modeled as a two-piece power law: /„ oc v^" 1 and /„ oc u~"' 2 with 
a break E\,. In the fit, the low energy spectral index f3± = —2 is adopted, suggesting a self-absorbed 
prompt emission spectrum. fi% is adopted to be equal to the reported prompt 7-ray spectral index. 
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the existence of a supernova associated with GRB 090417B. 



5. Discussion 

All of the evidence presented in this paper points to GRB 090417B being a dark burst because 
of localized extinction in the host galaxy. The host galaxy has a redshift of z = 0.345, and is clearly 
seen at optical and ultraviolet wavelengths. Therefore Lyman-a absorption can not account for the 
missing optical afterglow. The X-ray light curve and spectrum during the first ~ 70 ks obey the 
conventional closure relationships and are consistent with synchrotron radiation from a relativistic 
fireball expanding into an external medium (although the structure of that medium can not be 
determined with the available data). Optical afterglows have been detected for many GRBs with 
similar X-ray properties, so it unlikely that the lack of an optical detection is due to unusual physics 
during the afterglow phase of the burst. Optical observations of GRB 090417B started 378 s after 
the BAT trigger, and follow-up observations were made from several observatories for several days 
after the burst. These observations went deep enough to have detected an optical afterglow if one 
had been present with a luminosity similar to that of other optical afterglows. Therefore, the dark 
nature of GRB 0904 17B can not be explained by a lack of observations. 

The dust hypothesis, however, can explain this dark burst. The X-ray spectrum is consistent 
with an H I column density of Xh ~ 2 x 10 22 cm -2 , which corresponds to an extinction of A v ~ 
11 mag along the line of site in the host galaxy, assuming a Milky Way extinction law as explained 
in § 12.21 This value is in agreement with the minimum extinction, as explained in § 14. 1\ that 
is needed to obscure the expected bright optical afterglow (Ay > 2.5 mag). An extremely high 
extinction (Ay ~ 15-40 mag) is also predicted by the dust scattering model, which is consistent 
with the X-ray data. Such a large amount of dust along the line of sight to the GRB naturally 
explains the lack of any detection at ultraviolet, optical, or near-infrared wavelengths while allowing 



a conv entional X-ray afterglow to be observed. It is also consistent with the result of iPerley et al 



(|2009l ) who found large mean extinctions for the putative host galaxies of several GRBs. Our 
result, that the extinction along the line of sight in the host GRB 0904 17B is likely to be at least 
ten magnitudes in the UVOT v band, provides strong evidence that local dust concentrations in 
the host galaxy are responsible for at least some low-redshift dark bursts. 

Unlike many GRB host galaxies the dust in the host galaxy of GRB 090417B appears to follow 
a Milky Way extinction law. In addition the host appears to be an L* galaxy with a star- formation 
rate that exceed ~ 1 A4q yr _1 . In other words, the host appears to have dust properties, and a 
luminosity that are similar to those of the Milky Way, and a star-formation rate that is consistent 
with what is seen in other GRB host galaxies. 

GRB 090417B is one of the few dark GRBs where the association with a host galaxy is secure. 
In most cases associations between a GRB and a galaxy have been made based on the probability 
of the nearest observed galaxy lying as close to the centre of the X-ray error circle as it does. There 
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are a few exceptions, such as GRB 000210, GRB 050713A and XRF 050416A. The host ga l axy o f 
GRB 000210 was identified from the locations of its radio and X-ray afterglow. iPiro et al.l (J2002) 
found a probability of a chance alignment of the galaxy and the burst of P c h = 0.016 and conclude 
that if this is the host then the dark nature of GRB 000210 is likely due to dust either at the 
location of the progenit o r or along the line of sight. XRF 050416A occurred at a low redshift 



(z = 0.6535 ICenko et al.l (J2005I )). Although technical ly dark an optical afterglow was detecte d and 



was consistent with the standard fireball model (e.g. JCenko fc Foxll2005l : iHolland et al.ll2007l ). The 
presence of an optical afterglow and the low line-of-sight extinction suggests that this burst should 
not be considered to be a dark GRB. T he line-of-sight extinction for XRF 050416A is Ay < 1 mag 
([Holland et al.l 120071 : IPerlev et al.l 120091 ). which is consistent with what is seen for most optically 
detected GRB afterglo ws. IPerlev et all (0 009) note that GRB 05071 3 A also only barely qualifies as 
a dark burst under the iJakobsson et al.l (J2004J ) definition. The internal line-of-sig ht extinction for 
this burst is unknown. However, an optical afterglow was seen (jWren et alj|20 05), so it is unlikely 
that the extinction was large. The latter two of these bursts are techni cally dark based on th e 
spectral slope between the X-ray and optical regimes {Pox) criteria of IJakobsson et al.l (|2004l ). 
However, in these cases optical or infrared afterglows were detected, suggesting that these GRBs 
were borderline cases of dark bursts. 

Further evidence that dust is responsible for GRB 090417B being dark comes from the late-time 
softening of the X-ray spectrum. This softening is consistent with wh at is predicted from the dust 
scattering model (jShao fe Daill2007l ; IShao et al.ll2008l ; IShen et al.ll2009l ) , and requires approximately 
15-40 mag of extinction at optical wavelengths. The effects of dust scattering do not become 
apparent in the X-ray spectrum until ~ 10 5 s after the burst. Since X-ray spectra are usually 
extracted from early-time data (when there are a large number of photons) it is possible that the 
softening seen for GRB 090417B has been present in the late-time X-ray spectra of other dark 
bursts, but has not bee n seen due to low count rates. For example, the standard catalogue of 
Swift/XRT spectral fits (JEvans et al.ll2009l ) uses data taken within ~ 4.3 ks of the BAT trigger, 
well before a softening due to dust scattering is expected to appear. 



Per ley et al.1 (120091 ) point out that the host galaxies of dark GRBs are not unusually reddened 
relative to the hosts of optically bright GRBs. Their analysis assumes that their statistical as- 
sociation of bursts with the nearest detected galaxy on the sky reveals the true hosts of these 
bursts. It is uncertain if this is valid for all of the GRBs in their sample. However, it is likely 
that at least some (if not all) of these galaxies actually did host the associated GRBs. Although 
some of the hosts of dark GRBs have been f ound to be unusua lly dusty compared to the hosts 
of optically-bright GRBs (e.g., GRB 030115 (JLevan et al.ll2006l )). there is no evidence that the 
hosts of dark GRBs have systematically higher mean extinctions than the hosts of optically bright 
GRBs. However, the mean extinction of the host is probably not the factor that determines if a 
particular burst is dark or not since any extinction due to material that is not along the line of 
sight to the burst will not affect our observations of the afterglow. What is important in making a 
burst optically dark or bright is the extinction along the line of sight. This extinction can be due 
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to dust at the location of the GRB, although it has been suggested that the ultraviolet and X-ray 
emission from GRBs can destroy circumburst dust out to a few tens of pc from the progenitor 



(jWaxman &: Drainel 12000 : iFruchter et al.l l200ll ) . One way to get a dusty environment around a 
GRB progenitor is the stellar wind from a Wolf-Rayet progenitor. Collisions between wind-driven 
shells of material ejected during various stages of the Wolf-Rayet star's evolut i on can trigg er the 



formation of large amounts of dust (e.g. I Williams et al. 



1990 



Shao et al. 



2008; 



Cherchnefi 



20091 ). 



The radius of a Wolf-Rayet wind bubble is consistent with with the distance from the progenitor 
that we find for the dust along the line of sight to GRB 090417B (tens of pc). Alternately, the 
extinction could be due to unrelated dust in the host that just happens to lie along the line of 
sight, such as a dusty star-forming r egion or a giant mo lecular cloud. Giant molecular clouds have 
diameters of approximately 100 pc (JMizuno et al.ll200ll ). but often contain dense cores where star 



formation occurs. A chance alignment of one of these cores along the line of sight to GRB 090417B 
could account for the high extinction inferred from the X-ray data. 

For GRB 090417B the late-time evolution of the X-ray spectrum can be explained if there is 
a sheet of dust at a distance of 30 < Rd < 80 pc from the burst. This suggests that GRB 090417B 
may have been dark because of either a dusty e nvironment within 30-80 pc of the progenitor 



(e.g. IWaxman 8z Draind |2000| ; IFruchter et al.ll200ll ). or a chance alignment of the progenitor with 



a regi on of high extinction in the host galaxy The Galactic extinction model of iDrimmel et al. 



(J2003I ) suggests t hat roughly 25% of the sight lines through the Milky Way have Ay > 1 mag. 



The modelling of lUpdike et al.l ([2009) suggests that most sightline through a galaxy like the Milky 
Way, as seen from within the galaxy, have Ay < 1 mag. However, their work does not include 
dumpiness in the dust distribution, which would lead to a larger fraction of sightlines having higher 
extinctions. This suggests that a non-negligible (albeit uncertain) fraction of GRBs may be located 
along a sightline with more than approximately one magnitude of extinction in their hosts. When 
high-redshift dark bursts are accounted for the fraction of dark bursts is roughly comparable of the 
estimated fraction of highly extincted sightlines, which suggests that obscuration by dust along the 
line of sight may be responsible for some dark GRBs. 



6. Conclusions 



We have shown that the host of the dark GRB 090417B is very likely to be SDSS J135846.65+470104.5, 
an L* galaxy at z = 0.345. This galaxy has an overall star-formation rate that is at least as great 
as that of the Milky Way, and an overall F-band extinction of Ay ~ 3.5 mag, which is dustier than 
typical GRB host galaxies, but not greatly different from the Milky Way. 

X-ray observations of GRB 090417B show the normal temporal and spectral behaviour seen 
in the X-ray afterglows of many Swift GRBs. The data up to approximately 70 ks after the BAT 
trigger are consistent with a relativistic fireball expanding into either a homogeneous ISM or a 
wind-stratified external medium. The cooling frequency is below the X-ray band between 9.1 ks 
and 29.6 ks, and the /c-corrected isotropic energy is E[ so > 6.3 x 10 51 erg. The physics of this burst 
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do not appear to be unusual, and thus are unlikely to explain the dark nature of the burst. 

After approximately 70 ks the X-ray spectrum b ecomes signi fi cantly softer. We find that this 



can be explained using the dust scattering model of IShao fc Pail (|2007l ). The observed late-time 
spectral evolution can be produced by a sheet of dust approximately 30-80 pc from the burst along 
the line of sight. The model predicts an extinction of Ay ~ 15-40 mag along the line of sight. 
This is consistent with the H I column density measured from the X-ray spectrum. Therefore, we 
conclude that GRB 090417B is probably dark because of a dense, localized layer of dust along the 
line of sight between us and the afterglow. 

We acknowledge the use of public data from the Swift Data Archive. This paper is based 
in part on observations taken with the Nordic Optical Telescope, operated on the island of Santa 
Miguel de la Palma jointly by Denmark, Finland, Iceland, Norway, and Sweden in the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. The authors 
would like to thank the anonymous referee for their comments, which improved this paper. 
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